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Guided modes of multimode waveguides exchange power if the waveguide 
deviates in any way from its perfect geometry. The power exchange problem 
is studied for a multimode slab waveguide under the assumption that the 
power coupling is caused by irregularities of the core-cladding interfaces. 
The problem is treated by means of coupled power equations. The main 
result of this study is the realization that the power distribution versus mode 
number settles down to a steady state distribution if the waveguide is 
sufficiently long. The shape of the steady stale distribution depends on the 
correlation length of the function describing the core-cladding interface 
irregularities. For very short correlation length only the lowest-order mode 
carries an appreciable amount of power while the power carried by all 
the other modes is orders of magnitude smaller. For very long correlation 
length, on the other hand, all guided modes carry equal amounts of power. 
The steady state distribution is achieved regardless of the way in which 
the power was distributed over all the modes at the beginning of the guide. 
However, the total power in the steady stale mode distribution is dependent 
on the initial power distribution. 

I. INTRODUCTION 

Light communications systems using optical fibers as the guidance 
medium are presently being planned for two different modes of opera- 
tion. High-capacity systems are likely to be used with a laser as the 
light source and should be operated in the fundamental HE,, mode in 
order to minimize delay distortion that accompanies multimode opera- 
tion. For less ambitious, low-capacity systems excitation of the fiber 
with a light emitting diode appears more economical. However, the 
output of light emitting diodes cannot be used to excite a single fiber 
mode with high efficiency. A low-capacity fiber to be used with a light 
emitting diode must thus be designed to operate with many modes. 

Multimode optical fibers are not as easily characterized as single-mode 
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fibers. The power loss of such a fiber is usually not simply expressed 
by an exponential decay law but depends in a complicated way on the 
distribution of the power over the many modes. The present study is 
an attempt to describe the loss behavior of multimode optical fibers. 
We use the TE modes of the simplified model of a slab waveguide with 
the added requirement that there is no field variation or change in the 
slab geometry in the y direction of the coordinate system. This model 
makes it possible to describe the multimode waveguide rather simply. 
Even though it cannot directly be used to predict the loss behavior of 
round multimode optical fibers, it provides insight into the operating 
principles of multimode waveguides that can be used to obtain an 
understanding of the properties of multimode fibers of different shape. 
Our treatment of the multimode dielectric slab waveguide is based on 
coupled power equations. It has been shown in an earlier paper 1 that 
the coupled wave equations of a multimode optical waveguide 2 can be 
used to derive much simpler coupled power equations provided that 
the coupling mechanism can be described by a stationary random 
process with Gaussian correlation function. The coupled power equa- 
tions have the advantage that their coefficient matrix is constant, real, 
and symmetric. The system of coupled linear first-order differential 
equations can thus be solved by first finding eigensolutions with the 
common z dependence exp (-az). These can be used to express the 
general solution as a superposition of eigensolutions. This approach 
makes it clear that a steady state power distribution must exist. By 
allowing the field to travel far enough in the waveguide, so that all 
but the lowest-loss eigensolution has decayed to insignificant values, 
it is obvious that the distribution of power over the many modes 
assumes the shape of the lowest-order eigensolution regardless of the 
initial power distribution. The power loss of the steady state eigen- 
solution obeys a simple exponential law and can thus be characterized 
by a single number, the lowest-order eigenvalue of the eigensolutions 
of the power rate equations. 

The mechanism causing coupling between the many guided modes 
and of guided modes to the continuous spectrum of radiation modes will 
be assumed to consist of irregularities of the core-cladding interface. 
The coupling coefficients for this model have been evaluated in an 
earlier paper. 2 Any imperfection of the refractive index distribution 
and the slab geometry causes coupling between the modes. We choose 
the core-cladding interface irregularities because this coupling mecha- 
nism is of fundamental importance and because its properties are well 
understood. Mode coupling caused by irregularities of the refractive 
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index distribution will cause similar effects in some respects. However, 
there are differences that consist mostly in the dependence of the 
coupling process on the mode number of the coupled modes. 

The coupling coefficient for the core-cladding irregularities can be 
expressed as a product of a term that is independent of the length 
coordinate z but depends on the mode number times a z-dependent 
function that describes the actual shape of the core-cladding interface. 
This function /(z) is assumed to be a stationary random variable with 
a Gaussian correlation function that can be completely described by 
the rms deviation a of the core-cladding interface from a perfect plane 
and by the correlation length D. The same process that couples the 
guided modes among each other also causes each mode to lose power 
to the continuous spectrum of radiation modes. The interplay between 
coupling among the guided modes and power loss to radiation is respon- 
sible for the shape of the steady state distribution as well as for the loss 
associated with that stead}' state distribution. 

In order to spare readers not interested in the details of the theory 
the trouble of finding their way through the theoretical part of the 
paper, we present the results of the numerical analysis before the 
discussion of the details of the theory. 

II. RESULTS OF THE NUMERICAL ANALYSIS 

The theory has been evaluated for a slab waveguide with a core 
index of n, = 1.5 and a core-to-cladding-index ratio of n x /n 2 = 1.01. 
Most numerical results hold for a slab waveguide supporting ten modes 
corresponding to the value led = 82 (k = free-space propagation 
constant, d = slab half width). The only other case for which numerical 
values have been calculated corresponds to k<7 = 165 with twenty-one 
guided modes. It has been assumed throughout that the irregularities 
of the two core-cladding interfaces are statistically independent of each 
other but have the same rms deviation and the same correlation length. 

Figure 1 is a plot of the steady state distribution of the ten-mode slab 
waveguide. The steady state mode power is plotted versus mode number. 
Actually, only integer values of the mode number have physical mean- 
ing. In order to be able to display the mode power distributions for 
several values of the correlation length on one graph, the power values 
at the integer mode numbers were connected by straight lines. The 
label B { n u of the vertical axis refers to the lowest-order eigenvector of 
the eigenvalue problem [see equation (62) of the theoretical part]. 
These values are proportional to the power in each mode. They are 
normalized so that the squares of the power values for all ten integer 
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Fig. 1 — Steady state power distributions for different values of the correlation 
length D. The multimode guide carries 10 modes. 

mode numbers add up to unity. The most important aspect of Fig. 1 
is the shape of the steady state power distribution for different values 
of D/d, the ratio of correlation length to slab half width. For very long 
correlation length each mode carries an equal amount of power regard- 
less of the shape of the power versus mode distribution at the beginning 
of the guide. As D/d decreases more power is carried by the lower-order 
modes. For very small values of D/d (less than unity) essentially all 
the power is contained in the lowest-order mode. Figure 2 presents a 
similar graph for the case of twenty-one modes. The shape of the steady 
state distributions is essentially unchanged except that similarly shaped 
curves carry smaller D/d values showing that the number of modes 
does not affect the general behavior of the steady state distributions. 
The shape of the steady state distributions can be explained as 
follows. For long correlation length only the high-order guided modes 
lose power directly to the radiation field while the guided modes couple 
in such a way that only next neighbors exchange power. It is thus 
understandable that the power tends to equalize among all the modes. 
For very short correlation length all guided modes couple directly to 
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Fig. 2 — Same as Fig. 1 for the 2l-mode case. 



radiation. Higher-order modes lose power by this mechanism at a higher 
rate than lower-order modes. In addition, each guided mode couples 
to all the other guided modes. Since the lowest-order mode loses the 
least power to radiation it is the one that "survives" after all the other 
modes have lost nearly all of their power. In general, the correlation 
length of the random core-cladding interface irregularities cannot be 
chosen at will. However, for multimode operation one would hope for 
a long correlation length which makes it possible to transmit power in 
all the modes. Coupling with short correlation length forces the multi- 
mode fiber into single-mode steady state operation. 

Figure 3 shows the normalized steady state loss ad/(a\ 2 ) of the 
slab waveguide (the lines labeled i = 1) as functions of D/d. The lines 
labeled i = 2 represent the second eigenvalue of the eigenvalue problem. 
The important feature of Fig. 3 is the existence of a maximum as a 
function of D/d and the separation between the curves of the first 
(i = 1) and second (i — 2) eigenvalues. With the help of these two 
curves it is possible to estimate the region where steady state operation 
has been achieved. The loss parameters a lli and a <2) enter in the form 
exp ( — a l "z) as the first and second term of a series expansion [see 
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equation (62)]. When a i2) z > 4.6 we have exp (-a 2) z) < 10 2 so that 
the second term of the series expansion is becoming insignificant and 
steady state is essentially achieved. 

Figures 4, 5, and 6 show the way in which an initially uniform distri- 
bution of power settles down toward the steady state distribution. Three 
different values of correlation length were used. D/d = 0.01 is a suffi- 
ciently small value whose steady state distribution consists of only the 
lowest-order mode. The second mode carries only 10" 4 of the power of 
the first mode at z — * °o . The value D/d = 20 was chosen as an example 
for an intermediate correlation length. The steady state distribution 
in this case does not favor exclusively the lowest-order mode but assumes 
a shape in which higher-order modes carry decreasingly smaller amounts 
of power. The value of D/d = 35 is sufficiently large to produce an 
essentially uniform steady state distribution. 

The next three figures, Figs. 7, 8, and 9, show how the steady state 
distribution establishes itself if initially all the power is launched in the 
first mode. The last three figures, Figs. 10, 11, and 12, show similar 
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Fig. 3 — The normalized first two eigenvalues a (, '>d/(ff 2 k l ) (» = 1 and i = 2) are 
shown as functions of D/d for the 10- and 21-mode case. The lowest-order eigenvalue, 
t = 1, is the steady state power loss of the multimode waveguide. 
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Fig. 4 — Power distribution versus mode number for several values of normalized 
length along the guide for D/d = 0.01. 
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Fig. ~> — Same as Fig. 4, D/d = 20. 
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Fig. 6— Same as Fig. 4, D/d = 35. 
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Fig. 7 — Power distribution versus mode number for several values of the normal- 
ized length along the guide. Only mode 1 is excited at z = 0. D/d = 0.01. 
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Fig. 10— Power distribution versus mode number for several values of the normal- 
ized distance along the guide. Only mode 9 is excited at 2 = 0. D/d = 0.01. 




Fig. 11— Same as Fig. 10, D/d = 20. 
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Fig. 12— Same as Fig. 10, D/d = 35. 



plots for the case that all the power starts out in mode 9. The three 
values, D/d = 0.01, 20, and 35, have again been used. 

It may be of interest to know the ratio of the power remaining in the 
steady state if initially only mode 1 or mode 9 were excited. We call 
this ratio P(l, z)/P(9, z) and obtain P(l, z)/P(9, z) = 1.55 X 10 4 for 
D/d = 0.01, P(l, z)/P(9, z) = 9.75 for D/d = 20, and P(l, z)/P(9, z) = 
1.08 for D/d = 35. 

Tables I, II, and III show the ten eigenvalues of the steady state mode 
distributions together with the first eigenvector for the same three 
values of D/d. The lowest -order and the second eigenvalue appear also 
in Fig. 3. The other eigenvalues are given in the tables for the sake of 
completeness. It should be noted that the integer values in the left-hand 
column of these tables have different meaning for the eigenvalues and 
the eigenvector. The eigenvalues a ( " are ordered in increasing value 
and originate as the ten solutions of the eigenvalue problem [equation 
(60)] of a symmetric 10 by 10 matrix. The first eigenvector B\ l) belongs 
to the lowest eigenvalue a ( ". The subscript v is a mode label in this 
case. The eigenvector B[ x) is proportional to the steady state power 
distribution. Table I shows clearly that mode 1 carries the overwhelming 
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Table I — Eigenvalues and the First Eigenvector for 
D/d = 0.01 



i or v 


J- a M 

a 2 k* 


£„<»> 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


7.624 X 10- 6 
3.044 X 10-6 
6.845 X lO- 6 
1 .214 X 10"* 
1.891 X 10-*' 
2.711 X 10-* 
3.666 X 10"* 
4.743 X 10-* 
5.909 X 10"* 
7.049 X 10"* 


9.999 X 10-' 
8.486 X IO- 5 
7.160 X lO- 6 
6.789 X lO" 6 
6.630 X lO- 6 
6.546 X lO- 6 
6.497 X lO" 5 
6.466 X lO" 5 
6.444 X lO- 5 
6.424 X lO" 5 



Table II— Eigenvalues and the First Eigenvector for D/d = 20 



7 or v 


d (.) 
ff 2 k s 


B,<" 


1 
2 
3 

4 
5 
6 
7 
8 
9 
10 


5.175 X 10" 8 
1.402 X lO" 6 
4.576 X lO" 5 
1 .008 X 10-* 
1.951 X 10-* 
3.329 X 10-* 
5.217 X 10-* 
7.596 X 10-* 
1 .046 X 10- 3 
4.345 X 10" 3 


7.284 X 10-' 
3.799 X 10-' 
3.212 X 10-' 
2.770 X 10-' 
2.376 X 10-' 
2.011 X 10-' 
1.648 X 10-' 
1 .248 X 10-' 
7.475 X lO" 2 
3.129 X 10" 3 



Table III— Eigenvalues and the First Eigenvector for D/d = 35 



i or v 


<7*k 2 


fl F co 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 


6.415 X 10-" 
2.552 X 10~ 7 
1.112 X lO" 8 
4.251 X lO" 6 
7.789 X lO" 8 
1 .608 X lO" 6 
3.149 X lO" 5 
5.837 X lO" 5 
1 .036 X 10-* 
1 .765 X 10-* 


3.296 X 10-' 
3.294 X lO" 1 
3.293 X 10-' 
3.291 X lO-' 
3.290 X 10-i 
3.287 X 10-i 
3.279 X 10-i 
3.248 X 10-' 
3.094 X 10-' 
2.067 X 10"' 
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amount of power in the steady state if D/d = 0.01. The distribution of 
the first eigenvector for D/d = 20 appears also in Fig. 1. 

In order to obtain a feeling for the amount of irregularity of the 
core-cladding interface that will cause a certain amount of loss, we 
consider the case kd = 82 (10 modes). For X = 1 nm we then obtain 
d = 13 /xm for the slab half width. We now ask the question: What value 
of the rms deviation a causes 10 dB/km steady state radiation loss? 
The result is obtained from Fig. 3 or from Tables I, II, and III and is 
shown in Table IV. The tolerance requirements, arising from the need 

Table IV — Rms Deviation of Core-Cladding Interface 
Causing 10 dB/km Loss for kd = 82, X = 1 ^m 



D/d 


a/d 


a{nm) 


0.01 
0.3 
20.0 
35.0 


7.65 X 10-' 
3.26 X 10"' 
9.30 X 10"' 
2.63 X 10" 2 


9.94 X 10" 3 
4.25 X 10- 3 
1 .21 X 10-* 
3.42 X 10-' 



for keeping the steady state radiation losses low, are thus very stringent 
since the rms deviation of the core-cladding interface must be kept 
within fractions of micrometers. 



III. APPLICATION TO DELAY DISTORTION 

Multimode waveguides suffer from delay distortion that occurs 
because the modes contributing to the power transmission travel with 
different group velocities. Modes with a higher group velocity arrive 
at the receiver earlier than modes with a slower group velocity. A pulse, 
whose power is shared in some way by many modes, is thus distorted 
and lengthened by this effect. If the modes exchange power rapidly 
among each other this pulse lengthening effect of multimode waveguides 
can be substantially reduced. S. D. Personick 4 first pointed out the 
beneficial effect of tight mode coupling for the reduction of pulse delay 
distortion. For a two-mode waveguide Personick's results have been 
confirmed by a rigorous analysis by H. E. Rowe and D. T. Young. 5 
Our present work has some applications to the reduction of delay 
distortion by mode mixing. It is clear that if the coupling between the 
modes is strong, as would be desirable for delay distortion reduction, 

* A more rigorous discussion of pulse distortion in multimode waveguides will be 
published in a later issue of B. S.T.J. 3 
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the steady state power distribution is reached sooner. From Fig. 1 we 
see that only if the correlation length is large do many modes contribute 
to energy transport in the waveguide. As far as delay distortion is 
concerned it might appear advantageous to operate with a short corre- 
lation length forcing the multimode waveguide into essentially single- 
mode steady state operation. However, this method has the disadvantage 
that most of the power that is initially launched into higher-order 
modes is lost by radiation so that the waveguide suffers high transient 
losses. If a light emitting diode is to be used as the transmitter, single- 
mode operation is most undesirable. That leaves us only with the choice 
of a long correlation length (if indeed we have a choice) to reduce the 
power loss from high-order modes. In the limit of very long correlation 
length all the modes are excited equally strongly in the steady state 
distribution. If we can still provide strong coupling between the modes 
there is a chance that the power will be exchanged among all the modes 
making possible the reduction of delay distortion by mode mixing. 
Mode mixing takes place via coupling between nearest neighbors in 
case of long correlation length. The diffusion of power from mode 1 to 
the highest-order mode and vice versa is thus likely to be slow. 

We can get a rough idea of the "speed" with which the power travels 
from mode 1 to mode 9 (or from mode 9 to mode 1) from Figs. 9 and 12. 
It is apparent from both figures that it takes approximately z& 2 k 2 /d = 
10~" to 10" 7 before the mode at the other end of the mode spectrum has 
received an appreciable amount of power from the mode that is initially 
excited. The same power diffusion must, of course, take place for any 
other excitation of the modes. But the effect becomes observable when 
we launch all the power in one mode and watch how it redistributes 
itself over the other modes. This redistribution of power is part of the 
transient behavior that results in the steady state distribution. It is 
thus possible to estimate the distance that is required for one transit 
of power from mode 1 to mode 9 (or from mode 9 to mode 1) by looking 
at the second eigenvalue. We know that the steady state is reached as 
soon as the second term in the series expansion (of power in terms of 
steady states) becomes negligible compared to the leading first term 
[see equation (62)]. The second term is quite small when a 2) z = 2.3. 
We thus define a diffusion length L d by the relation 

L d = ff|r (1) 

L d is the distance along the waveguide that is required for the power in 
one of the modes at the end of the mode spectrum to transfer an appre- 



RADIATION' LOSSES IN MULTIMODE WAVEGUIDES 



443 



ciable amount of power (which is a somewhat undefined quantity) 
across to the mode at the other end of the mode spectrum. For D/d = 35 
we obtain from Table III and equation (1) L d d 2 \c/d = 9 X 10 6 . From 
Figs. 9 and 12 we see that this is indeed a reasonable estimate for the 
distance required for a power exchange between mode 1 and 9. The 
definition (1) allows us immediately to determine the steady state power 
loss that accompanies this power diffusion among the guided modes. 
The steady state power loss that occurs over a distance z = L d is 
given by 



» r — 9 q - 

Jj.l — 1.6 ,2, 



(2) 



Table V shows a number of values for a ll) L d for various correlation 
lengths. The penalty in radiation loss that must be paid for this mode 
mixing process is relatively high but it improves with increasing corre- 
lation length. Mode mixing via the next neighbor power exchange is not 
likely to be very effective in reducing delay distortion since only a small 
fraction of power traveling initially in one mode is transferred to the 
mode at the other end of the mode spectrum in the distance L d . One 
might expect that many such diffusion distances would have to fit into 
the overall length of the guide before delay distortion reduction by mode 
mixing becomes appreciable. Table V shows that a largo correlation 
length to slab half width ratio is required in order to keep the loss per 
distance L d small. Also shown in the table is the normalized exchange 
length L d . The numbers were computed for kd = 82, the ten-mode case. 
For delay distortion equalization it appears desirable to make L d much 
shorter than the total guide length L. If we choose L/L d = 100, for 
example, we compute from the last column of Table V for D/d = 40 
with L = 1 km, kd = 82, X = 1 /im, and d = 13 /mi the value a = 3.14 
jum for the required rms deviation of the core-cladding interface irregu- 
larities. This value is much larger than accidental irregularities need 



Table V — Loss Penalty a u L d and Normalized Power 
Exchange Length L d for kd = 82 



D/d 


a<"L,,(dBl 


5 2 k 2 . 


20 
30 
35 
40 


3.7 
0.33 
0.25 
0.064 


1 .64 X 1()» 
4.41 X 10 s 
0.02 X 10" 
3.01 X 10" 
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to be. It is thus conceivable that an optical fiber could be designed 
with an intentional core-cladding interface irregularity with long- 
correlation length for the purpose of reducing pulse delay distortion. 

IV. THE COUPLING COEFFICIENTS 

In Ref. 1 coupled power equations were derived from the coupled 
wave equations. The coupled wave equations have the form 

-P= tc»A/^'. (3) 

dz ^„ 

With the coupling coefficient written as 

c ril = KJ(z)- (4) 

and with the assumption that the correlation function of /(«) is Gaussian, 

</(*)/(* - u)> = *V"" D " (5) 

« > indicates an ensemble average), the coupled equations for the 
average power assume the form 

dP, = _ a ,p, + ^ a 2 D t I K» I 2 e' [in/2, "'^ )] \P. - P.). (6) 
dz m=i 

The term — a,P r was added to account for the radiation losses of the 
modes. Coupling coefficients describing the coupling between the guided 
modes of a slab waveguide caused by core-cladding interface irregular- 
ities were derived in Ref. 2. To obtain the coupling coefficient c v „ from 
our earlier work, we observe that equations (53) and (60) of Ref. 2 
correspond to a solution by perturbation theory of equation (3) for the 
special case that only the lowest-order even guided TE mode of the 
slab waveguide is excited. Comparison between the corresponding 
perturbation solution of (3) and equations (53) and (60) of Ref. 2 allows 
us to find 

c " Ba «[|PA|(i+T^)(l+^]* l/W ( } ()J " () 
The symbols appearing in (7) have the following meaning: 

d = core half thickness 

w, = index of refraction of core material 

n 2 = index of refraction of cladding material 

k = 2t/A = free-space propagation constant 

j3„ = propagation constant of mode v 
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<y, = (tf - ^k 2 ) J (8) 

k, = Kk 2 - #)* (9) 

a, = h K ' d for ^O- 2 ' 4 --- (10) 

(sin K r d for v = 1, 3, 5 

/(z) = distortion function of upper core-cladding interface (/(*) = 

indicates a perfect interface at x = d) 
h(z) = distortion function of lower core-cladding interface (h(z) = 

indicates a perfect interface at x = — d). 

The propagation constants of the even and odd guided TE modes are 
obtained with the help of (8) and (9) from the eigenvalue equations. 
We have for even modes 

tan*/* = — v = 0, 2, 4, ••■ (11) 

Km 

and for odd modes 

tan*,d = -— v = 1, 3, 5, ••• . (12) 

With the help of the eigenvalue equations, we can express (10) in the 
following form: 

H^ItT-V, f ° r ^°' 2 ' 4 - (13) 

( ( _ !)<>-»/*] (n, -n 2 )\ for .. 1|8> 5... . 

It is convenient to describe the guided modes in terms of a mode angle 0„ . 
We can introduce this angle by the equations 

k, — W|k sin 0„ , (14) 

/3„ = 7?,k cos 9 y . (15) 

Equations (14) and (15) represent the transverse and longitudinal 
components of the propagation vector of a plane wave in the core and 
are clearly compatible with (9). The guided mode can be represented 
as a superposition of two plane waves traveling inside of the core of 
the slab waveguide. ±0„ is the angle that these plane waves form with 
the waveguide axis. 

The function /(z) appearing in (4) is replaced with the sum of j(z) 
and h(z) in (7). Assuming that the two functions are uncorrelated, and 
assuming further that they have the same correlation function, we find 

([/(*) - {-\Y^h{z)]\j{z - u) - (-iy'"h(z - u)]) = 2crV ,u/D, \ (16) 
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By collecting all our results we can finally express the absolute square 
value of K y „ defined by (4) in the form 

. , 2 n 2 k 2 sin 2 0, sin 2 d„ , ir . 

2d V + ^A 1 + SJ cos '- cos ^ 

For small values of (nj/n 2 - 1) we have d, « 1. For modes far from 
cutoff we have, in addition, yd » 1. Under these conditions, (17) 
simplifies to the expression 

i K , r f . »mz. da) 

Far from cutoff we can approximate the solution of the eigenvalue 
equations (11) and (12) as follows: 

kA = (v+l)| for , = 0, 1,2,3,4 ••• . (19) 

The mode angles can then be expressed as 

• a 7r(" + 1) /on\ 

sin d y = (,-u) 

2n,kd 

Finally, it is important to know the largest mode angle that can occur. 
Mode guidance ceases to exist when the angle, that the plane-wave 
components of the guided mode form with the core-cladding interface, 
exceeds the total internal reflection angle C defined by 

?i! cos Q t = n 2 . (21) 

For (ni/n 2 — 1) « 1 we obtain approximately 



[<■ - s)] 



f.- 2U-? 



(22) 



Combining (20) and (22) allows us to find an approximate value for 
the number of modes N that the slab waveguide can support: 

2 - 1. (23) 



2 

N = - n t kd 

IT 



2I1-& 



V. RADIATION LOSSES 

In order to be able to evaluate the coupled power equations (6) we 
need convenient approximations for the radiation losses a, of the 
guided modes. 
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The radiation loss problem has been solved in a general way in Ref. 6. 
Equation (14) of Ref. 6 gives the radiation losses of the even and odd 
guide modes. 

a, = f' k (| F(fi, - 0) | 2 >/,(/3) dp (24) 

■'-n.k 

with 

I (q) = vh ~ "a )k n, sin 8, 



'h-'l cos 0,1 1 + -^ 
7„a 



[" p cos 2 ad _^ p sin 2 ad "1 , 

Lp 2 cos 2 o-rf + o- 2 sin 2 ad p sin 2 o-rf + a 2 cos 2 trrfJ 

and with the Fourier coefficient of the core-cladding interface function 
F(fi. - = -J= /" fW"*-' u dz. (26) 

Equations (10), (13), (14), and (15) have been used to express (25) in 
this form. In addition, (25) has been multiplied with a factor 2 to 
account for the fact that both core-cladding interfaces have irregularities 
(contrary to the assumption in Ref. 6) that are statistically independent 
of each other but have the same correlation function. There are two 
new parameters in (25): 

P = (n 2 k 2 - /3 2 )* (27) 

and 

a = (rc 2 k 2 - /3 2 ) ! . (28) 

The parameter (3 is the propagation constant (in z direction) of the 
radiation modes. Using (5) and assuming that L » D we obtain from (26) 

(| Ftf, - 13) | 2 ) = Vna 2 Dc- [,D/2U(, '- p)] \ (29) 

The loss expression (24) must be simplified before it can be used for 
our purposes. We are interested only in multimode waveguides with 
kd » 1. The functions sin ad and cos ad thus vary rapidly as functions 
of fi. It is impossible to obtain an approximation for all values of D/d. 
We begin by assuming D/d <<C 1. In this case, we can replace the exponen- 
tial function in (29) by unity and obtain 
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n 



] - n a 2 )k 3 n t sin 2 6, _ 2 D 



i\ u d 



2V^cos^(l+^) 

. r k [ _pMJ_*d + psin'ori "I ^ (3Q) 

./-n.k Lp 2 cos 2 o-d + a 2 sin 2 <rd p" sin 2 ad + «r cos" <rd-l 

Consider the terms in the integrand. The first term can be written 

p cos 2 ad _ p cos 2 ad ,^, 

p 2 cos 2 ad + a' sin 2 ad p" + (o- J — p 2 ) sin" o-d 

The sine and cosine functions oscillate rapidly while p is only a slowly 
varying function. The contribution of the second term in the denomi- 
nator is slight since this term vanishes when the cosine term in the 
numerator assumes its maximum. On the other hand, when the second 
term in the denominator is largest, the numerator is zero so that the 
value of the denominator does not matter. We can thus write to a crude 
approximation 

^ otfrf. (32 ) 

p 
The average value of the cosine square function is 1/2 so that we 
approximate further 

0, m i- (33) 

It appears that this approximation may be very poor at p = 0. However, 
by converting the integration variable from p to 0, we see that 

48 - -g dp (34) 

showing that there is no pole at p = 0. By an analogous argument we 
find that the second term in the integrand can also be approximated as 

ft«f. (35) 

2p 

The integral in (30) thus assumes the value 

f ,k (Gj + G 2 ) dp « f° - dp = r . (36) 

J- n ,k J-ntk P 

For D/d < (1/2/iakd) we obtain the following approximation for the 
radiation loss of the cth guided mode of the slab waveguide: 



RADIATION LOSSES IX MILTIMODE WAVEGUIDES 449 

Vir (n 2 - n 2 )w,k 3 sin 2 6, _ 2 D 

2 cos e ( l + ^) 

Next we try to obtain an approximation for large values of the corre- 
lation length, D/d » 1. The general expression for the radiation loss 
of the j/th guided mode follows from (24), (25), and (29): 

(n 2 - n^w.k'Vsin 2 B, , D 
a. = - — — a - 



2 V* 


('+£ 


cos 


5, 


f' k ,- 


■\(D/2)lfi.-fi)\ 


. 




'-„,., 


P cos 2 ad 







P cos 2 o-rf p sin 2 o-rf 1 

p" cos"' o-tf + o-"' sin" ad p"' sin 2 crrf + a 2 cos 2 o-rfJ 



(38) 



The exponential function under the integral sign decreases very rapidly 
with increasing values of jS, - for D/d » 1. Since the largest value 
that can assume is /3 = n a k only the immediate vicinity of the upper 
limit of the integration range contributes to the integral. In this region 
we have p « w 2 k. In order to be able to work out approximations for 
the case of large correlation length we must consider two more sub- 
divisions, the case that D/d is small enough so that the exponential 
factor under the integral sign in (38) varies slowly compared to the 
rapid oscillations of the sine and cosine functions and the opposite case 
where the exponential function decays to insignificant values within 
one cycle of the oscillations of the oscillatory functions. 

In the first case, slowly varying exponential function, we can consider 
P and a approximately constant over one cycle of oscillation except 
for the ad term appearing in the argument of the oscillatory functions 
and consider the average of the integral over one period 

J_ ('■ , .,' cos% f ... df 

— R. Ja. n" Pns" erf} -4- rr~ Bin nil 



/3 2 — p\ •'A' p" cos" ad + a sin ad 

cos 2 X 



2r J n 



dx « — — • (39) 



p" cos" X + a' sin X a + p 

And, similarly, for the second term of the integrand we find 



1 f ' p sin 2 ad 
~ZT / ^ • / , _l » i— , # « — 40 

Pi Jfi, p sill ad + a 



r_d_ 
-\- a' COS" 1 ad a + p 
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Since only small values of p can contribute to the integral, because of the 
rapid decay of the exponential function, we use the approximations 

1 1 

<t + p~ y/n\ — n 2 2 \a 

and 



(41) 



*- "* - fc (42) 

Using all these approximations and the change of integration variable 

#--! **-£*. < 43 > 

the integral in (38) can be approximated by the following expression: 



/ink 
>ntk 



'd/3 



-[(D/2) (0,-/3)1 ,_ 

e ap ~ , 2 , /-i 2 

i 2 



i 6 «*M ' — ' 2 /" 2 i 

a + p n 2 k VW] — w 2 



f 

«U(D/2)(0,-n,k)] 



du 



kVn? - W2 D 

= n 2 .^ , {l-erf|"f08,-n,k)]}. (44) 
kDVnl -nl I L.2 J J 

We have now finally obtained the result that the radiation loss of the 
j»th mode can be approximated by 

^ m ^kV^-njsin 2 !, .2^ _ erf [D % _ n>k) "J|. (45) 
d(l + Jj) cos $. 

The range of applicability of (45) is obtained by considering that we 
must require the exponential function in (44) to change only slightly 
over the range Ap corresponding to A(crd) = 2?r. This condition can be 
expressed as 

£^(^-n a k + ^r)«l. (46) 

2w 2 k V 2n 2 k/ 

The increment Ap is obtained from 

Ap-2r £ i (47) 

p a 
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The condition (46) must hold primarily for small values of p. We use, 
therefore, 

P = v Ap. (48) 

With 77 being a small number such as 2 or 5 and obtained from (46), 
(47), and (48) with a « (n? - t£)*k, 



(fi. - n£)d + {n] ~ n f kd 



< § < 1 t —^h — mi !*■ (49 > 

|^[C8, - n 2 k)d + W (W ' W2 nz) J(n? - ^) f J 

The left-hand side of this inequality follows from the requirement that 
the exponential function in (44) must drop to small values as p grows 
from to approximately (?ij — nl)*k. 

Finally, we obtain an approximation for very large values of D/d if 
we assume that the exponential factor in (38) decreases very appreciably 
over an interval corresponding to one oscillation period of the oscillatory 
functions. We can now use the approximation 

a = {n\ - nl)*k (50) 

treating a as independent of p. Expanding the factor that multiplies 
the exponential function in (38) in a power series in terms of p at p = 0, 
keeping only the first none-vanishing term, results in 



f 



£ -HD/2)(P,-P)l< 
(-1 6 



[ 



p cos 2 ad , p sin 2 ad 1 , 

p 2 cos J ad + a 2 sin 2 ad p~ sin 2 ad + a 2 cos 2 ad J 



' 2 2 — i ( c °t 2 <*d + tan 2 ad) 

n 2 (n l — n 2 )k 

fp 2 exp{-[f(^-n 2k + 2 -4)] 2 }.p 

'■' "■' k> ' (cot 2 ad + tan 2 arf) e - ,(D/2, " } '- n • k), ^ (51) 



In 3 



k\n;-n;)(p,-n 2 kyD 

If 0, = ?i 2 k, (51) becomes infinitely large. The approximation leading 
to the solution of the integral is violated in this case and the result 
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becomes meaningless. However, this violation of the applicability of 
the approximate solution of (51) can happen only for the highest-order 
mode and only if it happens to be directly at its cutoff frequency. Our 
approximation, if applied to this case, gives a loss value that is too large. 
Using too large a radiation loss for the highest-order mode affects the 
power distribution in all the other modes only slightly. The radiation 
loss of the highest-order mode is large in any case. Power coupled from 
the neighboring guided modes to this mode is lost rapidly. Using too 
large a loss value for this mode makes little difference to any of the 
other modes. We thus use (51) for all values of B, . 

A more serious violation of the applicability of (51) occurs if either 
the function cot ad or the function tan ad should become infinite or at 
least very large. In both of these cases the integral assumes the form 

/"° I. -l UV2M0r -W d p ^ -L f e -[<J>/a><*-n.k+<pV2n,k»l' ^ 

Jo B n 2 k Jo 



f ^ J JT -[<O/2H0,-n,k)l' 

~ D[n 2 k(B, - n 2 k)]* 
For very large D/d the radiation loss approximation is 

, 3 •_ 2 Q -[<D/2)(0,-n,k)P 

_ 2 n,k sin 0„e 



(52) 



2d[n 2 k(B, - n 2 k)]*(l + -^J cos 0, 

(cot 2 ad + tan 2 ad) for tan ad ^ 



kD 2 (/3„ - n 2 k) 

and cot ad ^ (53) 

(n 2 i — nl) for tan ad = 

or cot ad = 0. 

Equation (53) holds for values of D/d that are much larger than the 
D/d values in the range indicated in (49) . 

VI. THE EIGENVALUE PROBLEM 

Knowing the coupling coefficients and the radiation losses allows us 
to determine the power distribution in the multimode dielectric slab 
waveguide as a function of the distance z along the guide. Introducing 
the abbreviations 

h yu = VTa 2 D\K vll \ 2 e- [W2)(P '-^ ] ' (54) 
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and 

b, = £ K , (55) 

M = l 

we can write the coupled power equations (6) in the form 

^= -(«, + &,)P r + EMV (56) 

The trial solution 

P, = B,e~ az (57) 

converts (56) into an eigenvalue problem 

E [k„ - (a, + &„ - a)S ril ]B r = 0. (58) 

The coefficient matrix of this problem is real and symmetric as can be 
seen from (54) and the condition (11) of Ref. 1. This latter condition 
can be expressed in the form 

I #,„ P = I K,, | 2 . (59) 

The symmetry condition (59) follows also directly from (17). The 
eigenvalue oc is obtained from the eigenvalue equation 

I K, - (a, + 6„ - a) 5,„ | = 0. (60) 

The vertical lines in (60) indicate that the determinant of the matrix, 
whose v\i element appears explicitly, must be formed. The eigenvalue 
equation is an algebraic equation of order iV providing TV different 
solutions for the eigenvalues a (,) . The eigenvectors, whose elements 
are Bj°, are mutually orthogonal and will be assumed to be normalized, 

XJBW = 5„ . (61) 

The general solution of (56) can now be expressed as a linear super- 
position of the N eigensolutions, 

PM = £c,#V"\ (62) 

The expansion coefficients c, must be determined from the given power 
distribution at z = 0. With the help of (61) we obtain from (62) 

c, = EBj°J\(0). (63) 
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VII. CONCLUSIONS 

In this paper we have shown that coupling between the guided modes 
of a multimode waveguide causes the power versus mode number 
distribution to settle down to a steady state provided the signal is 
allowed to travel far enough in the waveguide. This steady state applies, 
of course, only to the CW case. For very long correlation length of the 
core-cladding interface irregularities the steady state distribution 
contains equal power in all the modes. For very short correlation length, 
on the other hand, only the lowest-order mode carries an appreciable 
amount of power, forcing the fiber into single-mode steady state opera- 
tion. 

The results of this paper have some application to delay distortion 
equalization. If the power carried by the guided modes is exchanged 
rapidly among them, the pulse distortion caused by the different group 
velocities of the modes is partially compensated. Coupling among the 
modes is of necessity accompanied by radiation losses. Effective pulse 
delay distortion equalization has a chance of working only if the corre- 
lation length of the core-cladding irregularities is long since the penalty 
paid in radiation loss becomes high for short correlation length. In 
addition, only for long correlation length do all the modes carry power 
in the steady state distribution. 

A detailed discussion of the numerical results and the properties of 
multimode waveguides is to be found at the beginning of the paper. 
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